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ABSTRACT 

Simple self-consistent models of galaxy groups and clusters are tested against the re- 
sults of high-resolution adiabatic gasdynamical simulations. We investigate two models 
based on the existence of a 'universal' dark matter density profile and two versions 
of the beta-model. The mass distribution of relaxed clusters can be fitted by phe- 
nomenological formulae proposed in the literature. Haloes that have experienced a 
recent merging event are systematically less concentrated and show steeper profiles 
than relaxed objects near the centre. The hot X-ray emitting gas is found to be in 
approximate hydrostatic equilibrium with the dark matter potential, and it is well 
described by a polytropic equation of state. Analytic formulae for the gas density 
and temperature can be derived from these premises. Though able to reproduce the 
X-ray surface brightness, the beta-model is shown to provide a poor description of 
our numerical clusters. We find strong evidence of a 'universal' temperature profile 
that decreases by a factor of 2 — 3 from the centre to the virial radius. We claim 
that the spherically-averaged profiles of all physical properties of galaxy groups and 
clusters can be fitted with only two free parameters. Numerical resolution and entropy 
conservation play a key role in the shapes of the simulated profiles at small radii. 
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1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
structures in the universe, and as such they have often been 
considered as a canonical data set for cosmological tests. As 
a first-order approximation for analytical studies, they can 
be described as spherically symmetric systems in which the 
baryonic gas is in hydrostatic equilibrium with an underlying 
cold dark matter (CDM) halo. 

Even for this basic model, the radial distributions of 
both types of matter are far from being well understood. 
During the last two decades, a great effort has been de- 
voted to investigate the mass distribution of CDM haloes 
by means of numeric al N-body simu l ations . The systematic 
study undertaken bv lNavarro et all lll997t hereafter NFW) 
concluded that the simulated density profiles can be fitted 
by a single two-parameter function, valid from galactic to 
cluster scales. 

Similar results have been found in several indepen- 
dent papers, although there is still some debate about 
the innermost value of the logarithmic slope of the den- 
sity profile and its dependence o n resolution. For instance, 
iFukushige fc Makind i 19971 1200 lD.lMoore et al] <1998Ul999t 
hereafter MQGSL) or lGhigna et alJ (Il998l l200ot find a cen- 



tral slope steeper t han t he NFW profi l e. On the other hand, 
IJing fc Sutol fcOOd) and lKlvpin et alJ (1200 ll) claim that the 
actual value might depend on h alo mass, me r ging history 
and substructure. More recently, iPower et alJ (§2003) point 
out that the logarithmic slope becomes increasingly shal- 
lower inwards, with little sign of approaching an asymptotic 
value at the resolved radii. 

The situation is even less clear for the baryonic com- 
ponent. Most gas in the intracluster medium (ICM) is in 
the form of a hot diffuse X-ray emitting plasma, where the 
cooling time (except in the innermost regions) is typically 
longer than the age of the universe. Adiabatic gasdynamical 
simulations have therefore been used to study the forma- 
tion and evolution of galaxy groups and clusters in d iffer- 
ent cosmologies fe.g.lNavarro et~al . 1995; Evrard et al Jl996t 
iBrvan fc NormarJl998tlEke et alll99Sft . This early work al- 
ready showed that the gas distribution is also similar in all 
objects, but more extended than that of the dark matter. 
Temperature profiles are found to decrease significantly with 
radius, but an isothermal core (or even increasing tempera- 
tures) of variable extent have often been reported. 

Recently, many numerical studies have focused on the 
effects of radiative cooling and non-gravitational heating on 
the final properties of groups and clusters of galaxies (see 
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e.g. iBialek et alJ l200ll: iDave et al.l 120021: iMuanwong et alJ 
120021: iBorgani et alJl200lt iTornatore et al J 12003^ " Entropy 
injection at high redshift and/or removal of low-entropy gas 
modify the radial distribution of gas at small radii, making 
the ICM even more extended than it would be in the purely 
gravitational case. Since cool systems are more sensitive to 
additional heating, non-adiabatic processes have often been 
invoked to explain th e discrepancy between the observed 
Lx — Tx relati on (e.g. [Ed ge fc Steward Il99lh and the self- 
similar scaling jKaiserl^8a? ~ as well as the entrop y excess 
detected in galaxy groups (e.g. iPonman et"aT1ll999l) . 

The extent to which additional physics can influence 
the cluster structure outsid e the cooling radius is unfortu- 
nately still unclear (se e e.g. iLewis et alJl200ol : iPearce et alJ 
l2000l:lMotl et alJl2003h . In this work, we are concerned with 
the theoretical predictions that can be made about the ra- 
dial structure of galaxy groups and clusters beyond self- 
similarity arguments, in an attempt to gain some under- 
standing of the purely gravitational case. We will consider 
that the ICM gas is in thermally-supported hydrostatic equi- 
librium with the CDM halo, which in turn dominates the 
mass of the system, 

1 dP(r) _ GM(r) ^ GA/ dm (r) 
p g {r) dr r 2 r 2 

where P = (p s kT) / (pm p ) denotes gas pressure, G is New- 
ton's constant, and M is the mass enclosed within radius r. 
p g and T are the gas density and temperature, k is Boltz- 
mann's constant, m p is the proton mass and p is the mean 
molecular weight of an ionised plasma of primordial compo- 
sition (p ~ 0.6). 

Several approximations are made in equation Q. The 
first one is spherical symmetry. Second, bulk and random 
motions within the gas are not included. And third, baryons 
can contribute a significant fraction of the total mass at large 
radii. These effects can lead to appreciable departures from 
equation Q in systems that are indeed in hydrostatic equi- 
librium. On the other hand, clusters of galaxies are evolving 
systems, and such a condition must not necessarily be ful- 
filled, particularly by dynamically young objects. 

Assuming a polytropic equation of state, 

P{r)<xpl(r) (2) 

a one-to-one correspondence bet ween mass, gas dens ity and 
temperature can be established dMakino et al.llT998l) . Note 
that isothermality is a particular case of a polytropic relation 
(i.e. 7 = 1). The physical origin of a polytropic equation of 
state is still a matter under investigation. A possible justifi- 
cation can be given in terms of the entrop y profile, related t o 
the mass accretion history of the system dVoit et alj [2003). 

A different approach to the radial structure of galaxy 
groups an d clusters is the the so-ca l led /3 -model, pio- 
neered by ICavaliere fc Fus co-Tejrrijmcjj[l976l) and widely 
used thereafter (see e'!g lRo^ati^t^dT]2002T for a recent re- 
view). This model assumes that the ICM gas follows a 'uni- 
versal' King-like density profile 

p g (r) =P0 [1 + (r/r c ) 2 ] (3) 

and that its temperature is approximately constant. Then, 
the observed X-ray surface brightness is given by the expres- 
sion 



S x (8)=S [l + (8/r c ) 2 ]- 3 ' +1/2 (4) 

where 9 is the projected radius from the observed X-ray 
centroid and So is the central surface brightness. 

The physical meaning of /3 is the ratio between gas 
temperature and the line-of sight velocity dispersion of 
the CDM component. The difference between this ratio 
(/9a P ec ~ 1) and the fits obtained by applying (3) to ob- 
served clusters of galaxies (/3fl t ~ 0.7) has been often re - 
ferred to as the '/^-discrepancy' (see e.g. ISarazinl Il986t) . 
None the less, it has long been know n both from observa- 
tions (e .g. |Nejrrnanii^^Arnaudl ll999F) and simulations (e.g. 
lBartelrriaJiri^^Stenim^t3^996T) that the /3-model does not 
provide an optimal description of the gas density, and that 
the best value of /Sat tends to increase with the outermost 
radius used in the fit. 

In this paper, we investigate the radial distribution of 
dark and baryonic matter by means of adiabatic gasdynami- 
cal simulations. We assess the validity of the assumptions of 
hydrostatic equilibrium and a polytropic equation of state, 
and compare the density and temperature profiles derived 
from these tenets with a sample of 15 numerical clusters of 
galaxies. 

The accuracy of the /3-model has been previously ad- 
dressed in several numerical stu dies (e.g. lEvrard et al.l(l996l : 
iBartelmann fc Steinmet3 Il996l) . focusing mainly on the 
mass estimates based on the observed X-ray surface bright- 
ness. In the present work, we consider four analytical cluster 
models. Two of them rely on the existence of a 'universal' 
CDM density profile, whereas the others assume a /3-model 
for the gas distribution. First, we attempt to obtain a self- 
consistent fit to our numerical data with each one of these 
models, noting that baryonic and dark matter profiles are 
not independent. Then, we try to estimate several cluster 
profiles from the simulated X-ray surface brightness. 

We would like to make special emphasis on the impor- 
tance of numerical resolution and accurate entropy conser- 
vation. We will show that the inner structure of the ICM de- 
pends critically on these issues. More specifically, we claim 
that they are responsible for the systematic differences ob- 
served between Eulerian and La grangian integration tech- 
niques (see e.g. lFrenk et alJ ll999). 

Section [5] describes our numerical experiments. Analyt- 
ical cluster models are summarised in Section |3] Results are 
compared in Section |1J and the impact on the analysis and 
interpretation of current observations is discussed in Sec- 
tion |S] We sum up our main conclusions in Section H2 



2 NUMERICAL EXPERIMENTS 

We have carried out a series of high-resolution gasdynamical 
simulations of cluster formation in a flat CDM universe with 
non-vanishing cosmological constant (f2 m = 0.3; £7a = 0.7; 
h = 0.7; a s = 0.9; fi b = 0.02 h~ 2 ). A thoro ugh description o f 
the numerical experiments can be found in lAscasibar! 12003 ) . 

Simulation s have been run with the parallel Tree-SPH 
code Gadget (Spr ingel et alJliboif) . We have used a new 
non-public version of the code, kindly provided by Volker 
Springel, i n which the entropy of SPH particles is explicitly 
conserved dSpringel fc Hernauistl2002h . Similar experiments 
have been accomplished with the Adaptive Refinement 
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Tree (ART) pure N-body code jKravtsov et alll997t) from 
the same initial conditions. In order to test the reliability of 
our results, one of the clusters ha s been re-simulated w ith 
the gasdynamical version of ART jKravtsov et al.ll2002l) . 

2.1 Cluster sample 

In a cubic volume of 80 ft -1 Mpc on a side, an unconstrained 
realization of the power spectrum of density fluctuations cor- 
responding to the ACDM model was generated for a total of 
1024 3 Fourier modes. The density field was then resampled 
to a grid of 128 particles, which were displaced from their 
Lagrangian positions according to the Zeldovich approxima- 
tion up to z = 49. Their evolution until the present epoch 
is traced by means of a pure N-body simulation with 128 3 
dark matter particles. 

Our sample comprises a total number of 15 clusters of 
galaxies, selected from this preliminary low- resolution ex- 
periment. Each object has been re-simulated with higher 
resolution by mea ns of the multiple mass technique (see 
iKlypin et aljl200ll . for details). First, we located the par- 
ticles in a spherical region around the centre of mass of the 
128 3 counterpart at z = 0. Mass resolution is then increased 
by using smaller masses in the Lagrangian volume depicted 
by these particles, including the additional small-scale waves 
from the ACDM power spectrum in the new initial condi- 
tions. 

We use 3 levels of mass refinement, reaching an effective 
resolution of 512 3 CDM particles (2.96 x 10 s h' 1 M Q ). Gas 
has been added in the highest resolved area only. The total 
number of particles (dark+SPH) in this area is greater than 
2x 10 6 for all clusters. The gravitational softening length was 
set to e = 2 — 5 h -1 kpc , depending on num ber of particles 
within the virial radius iPower et aUl2003ft . The minimum 
smoothing length for SPH was fixed to the same value as 
e. In total, 7 independent numerical experiments have been 
performed, running Gadget on a SGI Origin 3800 parallel 
computer at Ciemat (Spain), using 32 CPU simultaneously. 
The average computing time needed to run each simulation 
was ~ 8 CPU days (6 x 10 5 s). The same clusters have 
also been simulated with the N-body version of ART on the 
Hitachi SVR at the LRZ (Germany). 

Table Q displays the physical properties of our clusters 
at z = 0. Objects have been sorted (and named) according 
to their virial mass at the present day. Clusters J2 and K2 
are an exception to this rule, since they are two small groups 
falling into Ji and Ki, respectively. 

In order to study the effect of mergers and close encoun- 
ters, we used the Bound Density Maxima galaxy fin ding 
algorithm (see e.g. lKlypin et alll999l : ICoh'n et all 19991) . We 
label as major merger any cluster in which we are able to 
identify a companion structure inside -R200 whose mass is 
greater than 0.5 M200; if the most massive companion is in 
the range [0.1 — O.5M200], the object is classified as a mi- 
nor merger, otherwise, we consider it a relaxed system in 
virial equilibrium. The results of this classification scheme 
are quoted in the second column of Table Q Clusters named 
Ji and Ki are relatively close pairs, but they are separate 
enough (~ 2 — 3 Mpc) not to be considered as mergers. 

We find a remarkable amount of dynamical activity in 
our randomly selected sample. Only 6 objects have been cat- 
alogued as relaxed clusters, whereas 5 fall into the category 



Cluster 


State 


7?200 


M200 


-^200 


T200 


^200 


A 


Minor 


931 


18.96 


70.06 


2.000 


2.873 


B 


Minor 


871 


15.57 


20.30 


1.860 


2.620 


C 


Relaxed 


871 


15.53 


42.65 


1.958 


2.810 


D 


Minor 


771 


10.77 


11.49 


1.301 


1.642 


E 


Relaxed 


719 


8.74 


12.53 


1.287 


1.866 


F 


Major 


661 


6.79 


12.47 


1.059 


1.195 


G 


Major 


638 


6.10 


4.50 


0.844 


0.818 


H 


Relaxed 


618 


5.56 


16.60 


1.107 


1.614 


I 


Major 


581 


4.60 


2.34 


0.666 


0.666 


Jl 


Relaxed 


584 


4.67 


9.99 


0.937 


1.498 


Ki 


Relaxed 


557 


4.06 


4.98 


0.690 


1.043 


L 


Minor 


547 


3.84 


1.21 


0.624 


0.779 


M 


Major 


503 


2.99 


0.64 


0.545 


0.610 


K 2 


Minor 


497 


2.89 


3.08 


0.470 


0.808 


J 2 


Relaxed 


491 


2.77 


8.22 


0.673 


0.979 



Table 1. Physical properties ol the clusters at 2 = 0, defined 
at 200 times the critical density. Cluster name, dynamical state, 
72200 i n h~ 1 kpc, enclosed mass in 10 13 fo — 1 Mq, bolometric X-ray 
luminosity in 10 25 /i erg s — 1 , and average temperatures (mass and 
emission-weighted) in keV. Individual images of our 15 clusters 
can be found in Appendix 1X1 



of minor mergers and 4 have be en classified as major m erg- 
ing systems. As pointed out bv iGottlober et ail <l200lli . we 
find that the typical merging rate in groups is much higher 
than in the more massive clusters. 

Columns 3-7 display the bulk properties of each ob- 
ject. The subscript '200' refers to the overdensity with re- 
spect to the critical value, p c ~ 2.8 x 10 11 h 2 M Mpc -3 . 
Although the density contrast predicted in the spherical 
collapse model is closer to 100 for our ACDM cosmology, 
we chose this value for consistency with most observational 
studies. Since most of the mass is concentrated in the cen- 
tral regions, the difference between the two prescriptions is 
not large. 

Throughout this paper, we assume that X-ray emission 
arises from bremsstrahlung radiation only. The total lumi- 
nosity of a set of N gas particles is therefore computed as 



1=1 



where A = 1.68 x 10 17 erg s" 1 M, 



" 2 Mpc 3 



keV 



-1/2 



(•») 



For equal 



mass particles, the emission-weighted temperature is given 
by 



.3/2 



i = l P'^ 
V-^iV ml/3 

E<=i P iT i 



(6) 



These formulae have been used to compute the correspond- 
ing entries in TableQ as well as the analytical and numerical 
profiles of X-ray related quantities. 

2.2 Comparison with other numerical techniques 

Since the late 1980s a variety of techniques have been de- 
veloped to simulate gasdynamics in a cosmological con- 
text. In part inspired by the success of the N-body scheme, 
the first gasdynamical techniques were based on a parti- 
cle representation of Lagrangian gas elements usin g the 
smoothed particle hydrodynamics (SPH) technique (ILucvl 
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Il977l: iGineold fc Monaghanlll977l) . Soon thereafter fixed- 
mesh Eule rian methods were adapted llCen et al.l Il99dt 
ICenl Il992l i a nd, more recently, Eu lerian methods with 
sub- meshing (Brvan & Norman 1995), dcformable moving 
meshes l|Gnedinlll995l |Penlll995l ll99Sl) or adaptive mesh 
refinement CAMR. iBrvan et alJll995HKravtsov et aljEool) 
have been developed, as well as extensions of the SPH 
technique JShapiro et alJll996t Isiaringel fc Hernauisdl20o3 
ISerna et alJl2003ft . 

The Santa Barbara cluster comparison project 
jFrenk et al.lll99ST) attempted to assess the extent to which 
existing modelling techniques gave consistent results in 
a realistic astrophysical application. The formation of a 
galaxy cluster in a SCDM universe was simulated with 
12 independent codes, seven of them based on the SPH 
scheme and five on Eulerian methods. The properties of 
the CDM component were encouragingly similar, most 
discrepancies arising from small differences in timing (there 
was a merging event at z ~ 0). Less agreement was found 
in the gas-related quantities, but usually still within 10— 20 
per cent. Only the X-ray luminosity showed a strong 
dependence on the resolution of the different codes, with a 
spread as high as a factor of 10. 

One of the most remarkable findings was a systematic 
trend in the temperature profiles obtained for the inner re- 
gions (r 100 kpc). Near the centre, SPH codes generate a 
flat (or even slightly declining inwards) temperature profile, 
while grid codes produce temperature profiles that are still 
rising at the resolution limit. The entropy profile in SPH 
codes decreases continously towards the centre, while grid 
codes develop an isentropic core at sm all radii. 

Several authors (e.g . iHernauisd 

Nelson fc Papaloizoul Il99l 11994 : ISerna et alJ Il99rj. 
Springel fc Hernouistl 120021) have pointed out that 



1993; 



2003; 



important shortcoming of conventional SPH formulations 
is the poor conservation of entropy when a low number 
of particles is used. We have run several test simulations 
of the Santa Barbara cluster to study the effects of nu- 
merical resolution and entropy conservation. The st andard 
implementation of Gadget iSpringel et al.l 1200 it) gives 
simila r results as the SPH-based codes in iFrenk et alJ 
(1999), while the temperat ure and density profiles of the 
entropy-conserving version (ISpringel fc Hernquisdl2002D arc 
closer to the results of the Eulerian c ode by IBrvan et alJ 
lll995t) up to the resolution limit (see lAscasibaJ l2003t for 
details). The innermost radius r m i n was defined according 
to the most restrictive of the following criteria: 

(i) 200 dark matter p articles jKlypin et al1l200 ill . 

(ii) 100 gas particles feorgani etal*ll2002l) . 

(iii ) 3 times the gravitational softening iPower et alJ 

Usually, the resolution limit is set by the last two con- 
ditions. Excessively low values of the smoothing length can 
lead to the formation of cold compact groups of SPH parti- 
cles that decouple from the surrounding medium. Although 
Eulerian simulations are not affected by this problem, an 
unphysical temperature drop is often found in the cores 
of clusters simulated with SPH-based codes at low and 
medium resolutions fe.g. [ Eke et al.lll998t IFrenk et"al]ll999t 
iMathiesen fc Evrardll200ll) . When the smoothing length is 
properly set, a high number of particles is required to notice 
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Figure 1. Simulations of Cluster A performed with the standard 
SPH (dashed lines) and entropy-conserving (solid lines) imple- 
mentations of Gadget, compared to the results of the Eulerian 
code ART (crosses). Top left: Dark matter density. Bottom left: 
Gas density. Top right: Gas temperature. Bottom right: Gas en- 
tropy. 



the effects of numerical entropy loss, because the entropy 
profile artificially flattens at r ^ r m j n due to the soft ened 
density and temperature (see e.g. iBorgani et aT] |2002). 

The maximum resolution allowed by the Santa Barbara 
initial conditions is only 256 3 particles, and the compar- 
ison between different runs is further complicated by the 
merger at z ~ 0. Therefore, we decided to test the gasdy- 
namical integration scheme with the most massive of our 
clu sters. The results of the Eulerian code ART for Cluster 
A jNagai fc Kraytsovll2002h are compared with both imple- 
mentations of Gadget in Figure 

As expected, CDM density profiles agree within a few 
percent. Moreover, we find that the mass distribution in 
the pure N-body ART runs is almost identical (AM/M ~ 
10 per cent) to that obtained in Gadget, once rescaled to 
account for the slightly different values of Sldm due to the 
presence of baryons. Only clusters F, I, J and K display 
bigger discrepancies (~ 40 per cent) related to offsets in 
timing (note that all these systems are major mergers). 

Differences between both implementations of Gadget 
are mostly evident in the gas distribution at small radii 
(~ 20 kpc). The density profile is steeper in the standard 
formulation of SPH, whereas the temperature is systemat- 
ically higher when entropy conservation is enforced. There 
is an excellent agreement between the gas density found in 
ART and the entropy-conserving version of Gadget, but 
we find a lower temperature in the latter throughout most 
of the cluster. 

The shape of the temperature profile, though, is simi- 
lar in the two codes, and noticeably different from the ap- 
proximately isothermal structure found in the standard SPH 
version. This affects the gas entropy near the centre, which 
shows a 'floor' value of 30-40 keV cm" " 2 at r ~ 20 h' 1 kpc 
in both ART and entropy-conserving Gadget runs, in con- 



The radial structure of galaxy groups and clusters 5 



trast with most SPH implementations (including the pub- 
lic version of Gadget), which tend to produce a decreas- 
ing entropy profi le down to the resolution limit (see e.g. 
iFrenk et alJll99Sh . 

The agreement between the Euler ian code ART and the 
explic it entropy-conserving scheme bv lSpringel fc Hernauisd 
(2002) provides encouraging support to this new SPH for- 
mulation. Furthermore, it constitutes a sound piece of evi- 
dence that the conventional scheme suffers from severe en- 
tropy losses, with important consequences on the shape of 
the profiles in the central regions of groups and clusters of 
galaxies. We conclude that a careful modelisation of adia- 
batic processes is an essential requisite prior to the inclusion 
of additional physics, such as cooling or preheating of the 
intergalactic medium at high redshift. 



3 ANALYTICAL MODELS 

We compare our numerical data with self-consistent analyt- 
ical cluster models based on the assumptions of hydrostatic 
equilibrium and a polytropic equation of state. As a first 
step, we check whether these approximations are consistent 
with our simulated sample of galaxy groups and clusters. 

Once the reliability of the basic assumptions is tested, 
we focus on two different kinds of models: one is based on 
the existence of a 'universal' dark matter density profile (e.g. 
NFW or MQGSL) and the other relies on a 'universal' gas 
density (i.e. /3-model). In total, we test four different descrip- 
tions of the radial structure of galaxy groups and clusters. 
Hereafter, we will refer to them as NFW, MQGSL, /3-model 
(BM), and polytropic /3-model (PBM). 



3.1 Basic assumptions 

The extent to which equations Q and @ hold for our nu- 
merical clusters is measured in Figure [5] On the top panel, 
we plot the ratio between the gravitational term GM/r 2 
and the pressure gradient p~ 1 dP/dr. For a gas in thermally- 
supported hydrostatic equilibrium (i.e. neglecting infall, an- 
gular momentum and turbulence), both quantities should 
be equal to satisfy equation Q. Clusters that have been 
classified as relaxed on dynamical grounds can be consid- 
ered in hydrostatic equilibrium with an accuracy better than 
10 per cent. For minor mergers, this assumption holds only 
marginally (~ 25 per cent). Clusters undergoing a major 
merger display an even larger scatter, particularly at small 
radii. 

Deviations from hydrostatic equilibrium arise mostly 
from the contribution of kinetic energy to the gravitational 
support of the ICM gas. In relaxed clusters, thermal en- 
ergy dominates the total energy budget. Velocity dispersion 
increases substantially during merging events, and the gas 
needs some time to dissipate this extra kinetic energy into 
heat. Since relaxed systems constitute only 40 per cent of 
our sample, some caution must be kept in mind whenever 
equation Q is applied to a real cluster ensemble. We note 
as well that infall and random motions are important in 
relaxed clusters for r ^ 0.8i?2oo- 

Concerning the equation of state, the most straightfor- 
ward way to check whether the ICM follows a polytropic 
relation is to compute the effective value of the polytropic 
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Figure 2. Hydrostatic equilibrium and polytropic equation of 
state in our numerical sample. Objects have been classified ac- 
cording to their dynamical state as explained in l2.ll Solid squares 
correspond to the average over all systems in each category, and 
error bars are used to indicate one-sigma deviation of individual 
profiles (dotted lines) around the average. Top panel: Test of the 
hydrostatic equilibrium assumption. Bottom panel: Polytropic in- 
dex as a function of r. 




Gas density [h 2 M. Mpc 3 ] 

Figure 3. Histogram of individual gas particles of Cluster A 
in the p g — T plane. Contour lines are drawn at 10, 100 and 
300 particles per bin. Dashed line shows the expected slope for 
7 = 1.18. 



index from the spherically-averaged density and tempera- 
ture profiles, 

1[T) 1+ dlog[p g (r)] {l) 

This quantity is shown on the bottom panel of Figure [5] 
Albeit some scatter, gas in our relaxed clusters and minor 
mergers is consistent with a polytropic equation of state 
with 7 ~ 1.18 (dashed line). An isothermal profile (7 = 1) 
can be confidently ruled out from these data. A density- 
temperature histogram of the gas particles of Cluster A is 
plotted in Figure [3] We see that a polytropic equation of 
state reflects not only the mean behaviour of the ICM gas, 
but also describes the actual ratio between the density and 
temperature of individual mass elements. Although only one 
of our clusters is shown in Figure El the plot looks similar 
for any other object. 

Gas in small haloes is typically much colder than the 
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intracluster medium. Cold gas clumps corresponding to in- 
falling galaxies can be distinctly appreciated in Figure^] but 
their mass is too low to alter the effective polytropic index of 
Cluster A. However, the profiles become very noisy beyond 
the virial radius, as well as in major mergers. A constant 7 
is not an accurate approximation in these cases, where the 
amount of substructure within a spherical shell is large. 

3.2 Cluster models based on a 'universal' CDM 
density profile 

As po inted out in iMakino et alJ il998l) and ISuto et "ail 

( f998), the gas density and temperature profiles can be de- 
rived analytically from the CDM mass distribution under 
the assumptions of hydrostatic equilibrium and a polytropic 
equation of state. In the present work, we will consider the 
functional forms proposed by NFW, 



JNFW / \ 

Pdm (r) ~ {r/r B ){l + r/r s y 
and MQGSL, 

„MQGSL (r) 



(8) 



(9) 



(r/r„)3/ a [1 + (r/r m )3/2] 

Gas density and temperature profiles can be computed 
by substituting these mass distributions into equation Q 
a nd combining it w ith equation J5J|. Following the notation 
of ISuto et alJ ( 19981) . the gas temperature is given by 



T(r) 
To 



l-Bi/i(r/n) 



(10) 



where the subscript T can be either V or 'm' for NFW and 
MQGSL models, 



/.(*) = 1 - 
and 1 



ln(f + x) 



tv 2 I 2x 1/2 

— = H = arctan 

3^ V3 



1 fx + 2x 1/2 
+ 3 n \ x - x 1 /' 2 + 1 



(11) 



V3 

1\ 2\n{l + x 



3/2^ 



3.r: 



In both cases, 

47rG , ^m p (7 - 



B 



l)pi^ 2 



ykT 



(12) 



(13) 



relates the central temperature To to the underlying dark 
matter distribution. The gas density can be easily obtained 
from the polytropic relation 



Ps( r ) 



T(r) 
To 



1/(7-1) 



(14) 



According to this simple picture, clusters of galaxies can 
be described as a function of 5 parameters. A characteristic 
density and radius define the properties of the CDM halo, 
whereas for the ICM we also need to specify the gas density 
and temperature at r — 0, as well as the polytropic index 7. 

However, equation 11311 relates the values of these vari- 
ables. The number of free parameters can be further re- 
duced by imposing additional constraints. For instance, 

1 Note the typos in formulas (15) and (16) of ISuto et all il99d) . 



iKomatsu fc Seliakl i200 ll) choose To and 7 by enforcing ap- 
proximately constant baryon fraction at large radii. In our 
case, we just impose vanishing density at infinity, i.e. B s — 1 
and B m — 3v3/(47r), and we compute the central gas den- 
sity from the condition that the baryon fraction never ex- 
ceeds the cosmic value, 



Pg(^max) 
Pdm(^max) ^dm 



(15) 



where r max is the radius at which the predicted baryon frac- 
tion reaches its maximum (2.71r s or 1.34r m ). With this pre- 
scription, 



NFW -1 r-i "b MQGSL 

Po -l- 5 !?^ — P s > Po 

* 'dm 



2;<,08-^p 



T^Pm (16) 



In principle, there is no reason why the baryon frac- 
tion cannot exceed the cosmic value somewhere within the 
cluster. Yet, our numerical results (see Section 14.31 below) 
indicate that the gas to dark matter ratio increases mono- 
tonically (i.e. r max = 00) up to fib/^dm- Since the analytical 
estimate obtained for NFW and MQGSL models does not 
tend to an asymptotic value, we decided to set the normali- 
sation at the maximum. The value of po varies by ~ 25 per 
cent for any choice between r max and -R200, but inaccura- 
cies in the assumption of thermally-supported hydrostatic 
equilibrium might bias a prescription based on large radii. 

The last parameter in this family of models is the poly- 
tropic index 7. According to our results (see Figure|5J, we set 
it to 7 = 1.18. Alt hough this value is cons i stent with recent 
obser vations (e.g. iMarkevitch et all Il998t ISanderson et alJ 
2003), it would be interesting to investigate its physical ori- 
gin, as well as whether it can be related to other parameters 
such as halo mass or concentration. However, such a study 
requires a much larger number of objects, over a broader 
mass range, in order to be statistically significant. 



3.3 /3-models 

Most observational studies do not explicitly rely on the as- 
sumption of a given CDM density profile. Instead, they fit 
the observed X-ray surface brightness with equation J3J , and 
then compute the dark and baryonic profiles by applying 
hydrostatic equilibrium and a polytropic (often isothermal) 
equation of state. 

To a great extent, the so-called '/3-discrepancy' is due 
to the fact that expression does not provide a good fit 
to t he gas distribution throughout the whole cluster (se e 
e.g. iNa varro et alJ 11995} [ Bartelmann fc Steinmetzl Il996l) . 
iMathiesen fc Evrardl (l200lD argue that line emission from 
cold gas in the outskirts of the cluster is entirely responsible 
for the observed value /3 ~ 2/3. 

Polytropic /3-models have also been proposed to account 
for the presence of a temperature gradient, but equation 
is still used to model the gas density (e.g. lEttorll200Cl) . In 
the general polytropic case, the surface brightness 



Sx(6>) = So 1 + (6/r c y 



(17) 



is similar to the conventional form Q for any reasonable 
value of the polytropic index (1 < 7 < 5/3). Therfore, the 
gas density inferred from a fit to the observed X-ray profile 
has virtually the same shape in both the isothermal (BM) 
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Model S /(ApgT 1/2 n ) (n bPi )/(n dln p ) T /(4*Gp,m pPi r?) 



NFW 
MQGSL 
BM 
PBM 



0.191 
0.0151 
1.571 
1.145 



1.51 
23.08 
1/3 
0.39 



0.153 
0.369 
6 

10.62 



Table 3. Relations between the characteristic radius r;, CDM 
density p;, gas density po, and gas temperature To. All quantities 
have been computed assuming 7 = 1.18, except in the BM (7 = 
1). For this model, j3 = 2/3; for PBM, /3 = 1. 



and polytropic (PBM) /3-models. The central surface bright- 
ness is given by 



So = ^fKp r c \T 



1/2 



rgi + - 1/2] 
r[3/3(i + V)] 



(18) 



Hydrostatic equilibrium leads to an underlying mass 
distribution 

3f3kT jr c 



M(r) 



-{r/r c y l + (r/r c y 



-3/3( 7 -l)/2-l 



(19) 



(20) 



G/im p 

which features a finite central density 

_ 9/?fcr o7 

Pdm(U) = p/s = —p, 2 

Therefore, the structure of the dark matter halo in the 
/3-model is related to gas temperature in a way analogous 
to the models described above. We can follow the same pro- 
cedure to normalise the central gas density, imposing that 
the maximum baryon fraction equals the cosmic value. We 
will consider a standard /3-model (7 = 1) with f3 = 2/3 and 
a polytropic form with 7 = 1.18 and (3=1. For these val- 
ues, the maximum baryon fraction occurs at r^^ x — > 00 and 



BM 

Po 



t 4.18r c , respectively, and hence 
1 fib 



3 fid 



-P/3 



PBM 

Po 



0.39-^ P/3 



(21) 



We quote in Table |5|the distributions of dark and bary- 
onic matter expected in our four analytical cluster models. 
Cumulative gas mass and projected quantities in NFW and 
MQGSL must be integrated numerically. Relations between 
the characteristic parameters are summarised in Table |H] 



4 RESULTS 

We have shown that our clusters are, to a fair extent, in 
hydrostatic equilibrium, and that a polytropic equation of 
state with 7 = 1.18 provides a good description of the hot 
diffuse component of the ICM. Under these conditions, mod- 
els of the radial structure of galaxy groups and clusters can 
be derived from phenomenological approaches to either the 
gas or dark matter density. In this Section, we compare these 
simple, self-consistent models with the spherically-averaged 
distributions of mass, gas density and temperature found in 
our sample of simulated galaxy clusters. 

For NFW and MQGSL, we fit the cumulative dark mass 
of each halo, while for the /3-models the gas mass is used in- 
stead. We compute the numerical profiles in 26 logarithmic 
bins between 0.05i?2oo and ifeoo- We set the lower cut-off 
well above our resolution limit (~ O.OI-R200) not only to 
avoid numerical effects but also to investigate how accurate 
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Figure 4. Best-fitting model parameters for our sample of nu- 
merical galaxy clusters (see text for details on the fitting pro- 
cedure). Solid squares represent relaxed objects, empty squares 
are used for minor mergers and stars for major merging systems. 
Contours indicate \Jx 2 l{dof) = 0.1. Dashed line in NFW panel 
plots c = 6.8TC,?,' . 



are the extrapolations of the analytical cluster models to- 
wards r = 0. 

The parameter space is explored by letting r2oo and c\ = 
7"20o/n vary uniformly in the range 0.25 < 7-200 < 1-25 h' 1 
Mpc and 0.1 < c ; < 16 (50 for BM). For the gas mass, 
the radius rio (enclosing a mean density of 10p c ) is used 
instead of r2oo- Both quantities are roughly equivalent for 
the baryon fraction assumed in this work. 

Results of the \ 2 minimisation are plotted in Figure |1| 
In most cases, the quality of the f it is only slightly better 
for NFW and MQGSL models {^Jx 2 /(dof) ~ 0.05 - 0.1) 
than for /3-models (0.05 — 0.15). Best-fitting values of r2oo 
are close to the actual R200 (quoted in Table obtained 
directly from the simulated mass profile. The values of no 
in BM (/3 = 2/3) are less well correlated with -R200 than 
those obtained for PBM (f3 = 1). 

Although the full sample is consistent with constant 
dark matter concentration over the restricted mass range 
probed, we notice a well-defined bimodal behaviour: while 
relaxed clusters follow the usual rela tion, approximated b y 
c = 186M2-Q0 1 for the NFW profile jBurkert fc Silklll999Fl . 
mergers are dramatically offset down, particularly low-mass 
systems. In the /3-models, we find that the best-fitting gas 
concentration even increases with cluster mass. In several 
systems, the core radius inferred by the BM is only slightly 
larger than the innermost radial bin included in the fit. 

An important fact is that dark and baryonic matter dis- 
tributions are expected to be strongly correlated, by virtue 
of the hydrostatic equilibrium equation and the constraint 
of a cosmic baryon fraction. All analytical models of galaxy 
groups and clusters studied in the present work have only 
two free parameters. According to our fitting procedure, gas 
density and temperature are genuine predictions (i.e. not 
fits) of the models based on a dark matter profile. Con- 
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NFW 



MQGSL 



BM / PBM 



PdmM x 1 (l + x) 2 

M dm (x) ln(l + x) - x(l + X)- 1 

p e (x) [x-Hnil + x)] 1 ^ 1 -^ 

T{x) aT 1 ln(l + a) 



x -3/2( 1 + :r 3/2 ) -i 

! '-» ' t 3 / 2 ) 

1/(7-1) 



1 



3 Mi 

3^ 



4ir 



/m(^) 



1 " ^/m(x) 



+ -B)x- 2 1 (l + x 2 )- 3S / 2 - 2 



i^a + x 2 )-^/ 2 - 1 



(l +a .3)-3^/2 
(1 + X 2 )" 3 - 8 / 2 



Table 2. Polytropic models of galaxy clusters in hydrostatic equilibrium, x denotes the radial coordinate in units of the characteristic 
radius of each model (r s in NFW, r m in MQGSL, and r c in /3- models). Dark matter density is scaled in terms of the characteristic 
density (p s , p m , and pp, respectively). Masses are in units of 4ir Pir?. Gas density and temperature are scaled by their central values po 
and Tq. The abbreviations fm(x), given by 1121 . and B = P(-y — 1) have been used. 
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Figure 5. Dark matter density profile of our numerical haloes. 
Black squares represent the average over all clusters, excluding 
major mergers. Error bars indicate one-sigma scatter of individual 
profiles. Left panel: CDM density, scaled by the best-fitting p; and 
rj of each object in the corresponding model. Analytical profiles 
are shown as solid lines. Right panel: Accuracy of each analytical 
profile. Vertical dotted lines mark the fitted region. 



versely, the same is true for the CDM density and gas tem- 
perature in the /3-models. 



4.1 Mass distribution 

The radial density profiles of our dark matter haloes are 
shown on the left panel of Figure |5] rescaled by their best- 
fitting characteristic densities and radii. We plot the ex- 
pected mass distribution in each model (see Table |5J as a 
solid line. On the right panel, we plot the accuracy of the 
analytical estimates, defined as 



Apdm _ Pest(f) 



,(r) 



(22) 



Pdm 

where p est and p nU m are the estimated and numerical pro- 
files, respectively. 

Not surprisingly, both NFW and MQGSL offer an ex- 
cellent fit to the dark matter distribution, but the core-like 




Figure 6. Logarithmic slope of the CDM mass (top) and density 
(bottom) profiles. Squares with error bars represent the average 
over relaxed clusters and minor mergers. Analytical cluster mod- 
els are shown as lines. 



predictions of the /3-models represent a very poor approxi- 
mation to the numerical profiles. The conventional /3-model 
has a much better accuracy throughout the fitted range, but 
the extrapolation to larger radii is more reliable in the poly- 
tropic version, thanks to its higher value of (3. 2 

Concerning the controversy about the inner slope of the 
density profile, it is apparent from Figure |^| that MQGSL 
formula is slightly more accurate than NFW near the cen- 
tre, albeit the scatter of individual profiles around the av- 
erage is somewhat higher. This issue is addressed in more 
detail in Figure H3 where we plot the logarithmic slopes of 
the CDM density and cumulative mass profiles, computed 
directly from the simulation data 

_ dlog[p d m(r)] dlog[M dm (r)] 



dlog(r) 



dlog(r) 



(23) 



In general terms, both NFW and MQGSL formulae are 
consistent with our results up to the resolution limit. When 
we divide our sample according to the dynamical state of 
each cluster, we find a well defined trend in the sense that re- 
laxed haloes display shallower central slopes (close to NFW) 
than merging systems (better described by MQGSL). As an 
extreme case, major mergers feature pure power-law mass 
distributions for more than one decade in radius, whose 
slopes are even steeper than the value proposed by MQGSL. 

2 In the isothermal case, though, dark matter density is insensi- 
tive to the actual value of this parameter (see Table |5J . 
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Figure 7. Same as Figure^ for the temperature profile. 



In agreement with IPower et alJ i2003Tl . we do not find any 
sign of an asymptotic slope for either relaxed clusters or 
minor mergers. More resolution is clearly required before 
achieving a firm conclusion on this matter. 

4.2 Gas temperature 

Mass determinations from X-ray emission in clusters usu- 
ally assume that relaxed clusters (i.e. morphologically sym- 
metric) are isothermal. However, this assumpti on has been 
questioned recently bv lMarkevitch et alJ a!998f) . who found 
evidence of a decreasing temperature profile for nearby clus- 
ters observed with ASCA. This result has b een confirmed 
for co ol clust e rs bv iFinogue nov et alJ <|200l|). On t he con- 
trary, IWhitd (|200(J> and llrwin fc Bregmanl i2000f> . using 
data from BeppoSAX and ASCA satellites, do not find any 
decrease of the temperature in a large collection of clusters. 
More recently, the analysis of BeppoSAX observations ac- 
complished bv lDe Grandi fc Molendil ((2002) concluded that 
temperature profiles of galaxy clusters can be described by 
an isothermal core followed by a rapid decrease. 

Numerical simulations of galaxy clusters support a sig- 
nificant decrease in the ICM te mperature from th e central 
regions to the virial radius (e.g. iFrenk et al.l|l999l) . Recent 
results from high resolution AMR gasdynamical simulations 
seem to indi cate that the tem perature profile has an uni- 
versal form jLoken et al J 12002). These authors proposed a 
simple formula to fit the projected emission-weighted tem- 
perature of their clusters: 

T p x (0) = T (l + 8/a x )- s (24) 

where Tp is the ce ntral temperature and a x a core radius. 
lLoken et alJ (120021) quote best-fitting values of these param- 
eters a x = r v ir/l-5 and 5 = 1.6. 

The average mass-weighted temperature profile of our 
sample of galaxy clusters (excluding major mergers) is 
shown on the left panel of Figure [7] The normalisation To 
has been computed analytically from the fit to the dark mat- 




B / R 200 

Figure 8. Projected emission-weighted temperature profile, 
scaled by T^ in Table HI ( black squares with error bars). Solid 
line is the best fit obtained with equation 124 1. dashe d line rep- 
resents the simu lation results of l Loken et al.l (2002). Observa- 
tional data from Grandi & Molcndi (2002) arc shown as tri- 
angles (filled for cooling flow and empty for non-cooling flow clus- 
ters). Dotted box enclos es the temperature profiles observed by 
iMarkevitch et al] ll998t) plus most of the error bars, solid box 
includes only the scatter of their best-fitting polytropic models. 



ter (NFW, MQGSL) or gas (BM, PBM) distribution of each 
cluster. As can be seen on the right panel, the theoretical 
prediction is accurate within ~ 15 per cent for all prescrip- 
tions, with the obvious exception of the isothermal /3-model. 
NFW and MQGSL models are virtually indistinguishable. 

In order to compa re with observation s, as well as with 
the numerical work of lLoken et al.l d2002tl . we plot in Fig- 
ure |H| the average projected X-ray temperature profile. We 
computed the emission-weighted temperature according to 
equation © in cylindrical shells of 3 ft -1 Mpc length, ori- 
ented along the main axes of the simulation box. The final 
profile of each cluster is given by the average over the three 
orthogonal projections, and the normalisation comes from 
the emission- weighted temperature, T^ , within -R200 (both 
quantities are quoted in Tabic 0. 

Our results a re in excellent agre e ment with the observa- 
tional estimate of IMarkevitch et alJ l|l99ftl) . who claim that 
the emission-weighted temperature is well represented by 
a polytropic /3-model (indeed, our data are almost identi- 
cal to their best-fitting profile). However, our profiles are 
only marginally consistent with the results reported by 
IDe Grandi fc M olcndi (2002). In particular, we do not find 
any evidence of the temperature flattening observed by these 
authors at r ~ 0.2i?200- 

The projected X-ray temperature of our clusters rises 
uninterruptedly until the innermost radial bin, in agreement 
with the results of Eulerian codes, such as ART (see Sec- 
tion ^, 21 or that employed bv lLoken et alJ ll2002l) . Equation 
$M fits well our numerical data, but we obtain rather differ- 
ent values of the core radius and the asymptotic exponent, 
a x = r2oo/4.5 and 5 = 0.7. Since both simulated tempera- 




ture profiles are consistent within the error bars, it seems 
that the best-fitting values of a x and S are rather sensitive 
to the details of the data. 

Never the less, we would like to stress that a decreasing 
temperature profile is a robust result of numerical simula- 
tions regardless of the integration scheme (either AMR or 
SPH), as long as entropy conservation is enforced in the 
latter. Furthermore, such behaviour is strongly supported 
by most recent observational measurements. The presence 
of an isothermal core in the observed galaxy clusters could 
reflect the need for additional non-adiabatic physics in the 
simulations. 

4.3 Gas density 

Finally, we compare the simulated gas density profile with 
the distributions expected in our four analytical cluster mod- 
els. It is evident from Figure |5] that the discrepancies are 
remarkably larger than for the dark matter density or ICM 
temperature profiles. 

It is somewhat surprising that the /3-models show the 
worst inconsistencies with the data, since they are specif- 
ically intended to fit the gas distribution. When /3 = 2/3, 
the asymptotic gas density falls too slowly, as r~ 2 . However, 
the central regions cannot be properly fitted if /3 = 1 is as- 
sumed. The final outcome is that the best-fitting value of 
P increases from one value to the other as larger radii are 
considered in the fit. 

Models based on a 'universal' dark matter profile pro- 
vide a fair description of the ICM density. None the less, we 
caution that in many cases the simulated gas density can dif- 
fer substantially from the analytical prediction, giving rise to 
the huge scatter (~ 40 per cent) observed in Figure|§] On one 
hand, the CDM distribution at small radii depends on the 
dynamical state of the halo. Near the virial radius, even re- 
laxed clusters deviate noticeably from thermally-supported 
hydrostatic equilibrium (see Figure |3J, and we also observe 
that the dark matter density is slightly higher than indi- 



cated by both NFW and MQGSL formulae (see Figure 0. 
More importantly, the polytropic index is not exactly the 
same for all clusters, which strongly alters the shape of the 
inferred gas density profile. 

Inaccuracies in the estimated dark matter and/or gas 
profiles lead to a significant mismatch between the analyt- 
ical and numerical baryon fractions. The local ratio of gas 
to CDM density is plotted in Figure fTOl as a function of ra- 
dius, rescaled by the cosmic value. None of the analytical 
cluster models is able to fit satisfactorily the simulated pro- 
file. The canonical /3-model is the only scenario in which 
the baryon fraction reaches an asymptotic value, but this 
is due to the equally wrong gas and dark matter densities. 
In the other models, the baryon fraction reaches a maxi- 
mum around three times the characteristic radius, but then 
it drops rather steeply. 

We used this maximum to find out the normalisation of 
the ICM density, po. However, the steep asymptotic slope of 
the gas density profile hints that more realistic models have 
to be considered in order to reproduce the radial structure 
near the virial radius. On the other hand, such a modelisa- 
tion (e.g. including departures from spherical symmetry, in- 
fall and turbulent motions) would introduce additional free 
parameters that would complicate the physical interpreta- 
tion of the results. 



5 OBSERVABLE CONSEQUENCES 

Analytical models of galaxy groups and clusters are a key 
ingredient in the interpretation of observational data, since 
they relate the X-ray emission of the ICM to the underly- 
ing gas and dark matter distributions. In particular, X-ray 
observations are usually restricted to the central regions, 
and these models allow the extrapolation of the spherically- 
averaged profiles up to the virial radius. In some occasions, 
they are also extrapolated inwards, in order to correct for 
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Figure 11. Best-fitting values of the central X-ray surface bright- 
ness, So, and the projected characteristic radius, Solid squares 
are used for relaxed systems, empty squares for minor mergers and 
stars for major mergers. Contours are drawn at yj \' 2 /(dof) = 0.2. 



the presence of a central cooling flow. Quite often, the virial 
radius itself is derived by assuming a particular model. 

Therefore, we would like to investigate in some detail 
to what extent the assumed radial structure can bias the 
conclusions drawn from observational studies. We construct 
X-ray surface brightness profiles in the same manner as we 
did for the projected emission-weighted temperature (i.e. as- 
suming pure bremsstrahlung radiation and averaging over 
the three axes). Then, we fit the simulated profiles with our 
four analytical models, varying the projected characteristic 
radius, 9i, and the central surface brightness, So. Results are 
shown in Figure ITTI 

Only data between 0.05 and 0.5 i?200 were used. The 
lower cut-off is similar to region excised in observational 
studies to avoid the cooling flow region. The upper scale is 
similar to that attained by X-ray data. Furthermore, we find 
that the emission at larger radii is significantly enhanced by 
small mass subhaloes, which manifest as peaks in the sur- 
face brightness. We discard these features by fitting only 
those bins that lea d to a monotonically decrea sing profile. 
As pointed out by Mat hiesen fc Evrardl i200ll) . the inclu- 
sion of emission lines would still increase the contribution of 
substructure to the total X-ray luminosity. 

The dependence of the fit on the choice of the outer ra- 
dius R out is plotted in Figure NFW and MQGSL models 
are more stable than /3-models, and their best-fitting 8i axe 
much closer to the actual characteristic radius n, computed 
from the mass distribution between 0.05 -R200 and -R200 (see 
Section|lJ|. The conventional /3-model offers a poor fit to our 
simulated surface brightness at large radii, due to the shal- 
low slope associated to j3 = 2/3. For several objects, the best 
fit is obtained with the lowest possible 9 C (i.e. a pure power 
law, Sx oc r~ 1,5 ). Since we assumed /3 = 1 for the PBM, 
a better agreement is found with the numerical results, al- 
though the best-fitting 8 C is systematically higher than the 
core radius inferred from the gas mass. 
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Figure 12. Dependence on the outer radius, -Rout, used for the 
fit to the X-ray surface brightness profile. Solid lines represent 
NFW model; dotted, MQGSL; short-dashed, BM; long-dashed, 
PBM. Left panel: Projected characteristic radius, divided by the 
three-dimensional estimation based on the CDM/gas distribution. 
Right panel: Goodness of the surface brightness fits. 

We assume that the characteristic radii of our clusters in 
each model are equal to their best-fitting projected equiv- 
alents, 9i. To obtain the characteristic CDM densities, as 
well as the central gas densities and temperatures, we must 
solve a system of equations, which in the /3-model is given 
by 1181 . 1201 and 1211 . For the models based on a 'univer- 
sal' dark matter density profile, the surface brightness is 
integrated numerically, yielding the central values quoted in 
Table |3] They must be combined with equations 1131 and 
1161 . substituting the appropriate value of Bi. 

A last interesting point to note is that we are not using 
any information about the X-ray temperature. In doing so, 
we could skip one of our constraints, such as the value of 
the polytropic index or the normalisation of the central gas 
density with respect to the CDM component. 

5.1 X-ray surface brightness 

The first question we would like to address is how well can 
our analytical cluster models fit the X-ray surface bright- 
ness. We show in Figure [T51 that the numerical profiles are 
fairly described by NFW or MQGSL, although the accuracy 
of the fits degrades dramatically at large radii. Never the 
less, we find that the difference between the analytical and 
simulated profiles is of the order of 30 — 40 per cent up to 
the virial radius. Some part of it is caused by the low gas 
density predicted by these models at large radii. The rest is 
due to the emission of substructure beyond 0.5 ifeoo- 

Consistently with observational results, the isothermal 
/3-model with (3 — 2/3 provides a fairly good fit to the X- 
ray surface brightness only for r ^ 0.3i?2oo- As was shown 
for the gas density, the outer parts are better approximated 
by f3 ~ 1. When we let f3 vary as a free parameter, the 
best fit to our data is obtained by (3 ~ 0.8. However, this is 
a 'compromise solution' that fails to properly describe the 
ICM at very small or large radii. 

5.2 Mass estimates 

One of the most important applications of hydrostatic equi- 
librium and a polytropic equation of state is the estimate 
of the cumulative mass profile of both dark and baryonic 
components, given the observed X-ray emission. We plot in 
Figure tTH the CDM mass distribution inferred for our cluster 
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Figure 13. X-ray surface brightness, averaged over all clusters 
except major mergers. Error bars indicate one-sigma scatter of 
individual profiles. Left panel: Surface brightness profile, scaled 
by the best-fitting So and 9; . Analytical cluster models are shown 
as solid lines. Right panel: Accuracy of each analytical profile. 
Vertical dotted lines mark the fitted region. 
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Figure 14. Dark matter mass profile, inferred from the X-ray 
surface brightness. Symbols as in Figure [T^l 



sample (excluding major mergers), scaled by the character- 
istic densities and radii derived from the fit to the surface 
brightness profile. 

The accuracy of the mass estimate is similar in all mod- 
els (~ 20 per cent errors), but the shape of the mass dis- 
tribution is substantially more accurate when a NFW or 
MQGSL 'universal' CDM density profile is assumed. The 
isothermal /3-model gives an acceptable estimate of the dark 
matter mass for 0.1 < r/i?200 < 1 because of the small val- 
ues inferred for the core radius. In the PBM, this quantity 
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Figure 15. Same as Figure H4l for the cumulative gas mass. 



is substantially higher and thus a 'core' in the dark matter 
distribution can be noticed for r ^ 0.3i?200- 

Cumulative gas mass estimates are plotted in Figure tlSl 
In this case, the small values of r c in the BM lead to 
extremely high central densities, whereas the PBM (with 
(3=1) predicts an excessively large core. Again, a profile 
closer to the numerical results can be obtained by treating 
j3 as a free parameter. The estimates based on a 'universal' 
CDM profile are only slightly more accurate. The scatter, as 
well as the overall shape, are greatly improved by consider- 
ing individual values of the polytropic index. 

Regarding the total baryon fraction, we can see in Fig- 
ure !16l that the /3-models can lead to severe overestimates of 
this quantity at small radii (r ^ 0.1 — 0.3 7?20o). Near the 
virial radius, though, uncertainties are of the order of 10 per 
cent. NFW and MQGSL models are accurate within ~ 20 
per cent for r ^ 0.2f?200 and ~ 50 percent as we move closer 
to the centre. 

Another (sometimes more important) source of uncer- 
tainty comes from the estimation of the virial radius itself. 
Furthermore, -R200 is not a fixed multiple of the character- 
istic radius in any analytical model of galaxy clusters, and 
hence the cumulative baryon fraction at a fixed overden- 
sity is extremely sensitive to halo concentration, particularly 
at low radii. Assuming a NFW model, the baryon fraction 
within -R200 can change by less than 10 per cent from c = 4 
to c = 9, but at 0.1 -R200, the ratio of gas to dark matter 
mass varies from .5 to .75 times the cosmic value. 



5.3 Entropy 

Entropy plays a fundamental role in clusters of galaxies be- 
cause convection acts as an entropy-sorting device, moving 
low entropy material to the cluster core and high entropy 
material to the outskirts. Gas density and temperature in 
hydrostatic and convective equilibrium are just manifesta- 
tions of the underlying entrop y distribution. The c laim of an 
entropy floor in small groups ( Po nman et al Il999l) has often 
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Figure 16. Same as Figure [T31 for the cumulative baryon fraction 
in units of the cosmic value. 



been regarded as evidence of preheati ng of the ICM at high 
redshift. Using a much larger sample. IPonman et alJ (120031) 
find that the observed profiles appear to be approximately 
self-similar apart from a normalisation constant. They ob- 
tain a best- fit behaviour S^O.lifeoo) oc T 0,65 for systems 
spanning a temperature range from 1 to 10 keV. 

The entropy p rofil e has recently been s tudied by 
iBorgani et alJ ll200lf) and lFinoeuenov et alJ (l2003f) by means 
of numerical simulations that include cooling and preheat- 
ing. With gravitational heating alone, these authors find 
a power-law entropy profile, consistent with other results 
based on the sta ndard implementation of SPH (see e.g. 
iFrenk et al]|l999h . As was discussed in Section [2.21 the sit- 
uation changes dramatically when entropy conservation is 
explicitly enforced, either in the SPH algorithm or in Eu- 
lerian schemes. Shock heating becomes substantially more 
efficient, leading to a much higher entropy near the cluster 
centre, as well as stripping low-entropy gas from the infalling 
galaxi es that otherwise w ould survive for several crossing 
times (IBorgani et al.ll200ll) . 

We plot our entropy profiles in Figure H71 As expected 
for a polytropic gas in hydrostatic equilibrium, they are not 
pure power laws for any reasonable choice of the 'universal' 
gas or dark matter distribution. Furthermore, the shape of 
the entropy profile does not depend systematically on the 
mass or temp erature of the object , in agreement with recent 
observations JPonman et alj |2003). NFW or MQGSL mod- 
els provide a much better estimate of the entropy distribu- 
tion than any version of the /3-model. However, the low gas 
densities predicted at large radii yield a slope at r ~ R200 
significantly steeper than the numerical data. 

Given our limited temperature coverage, the normali- 
sation of the entropy profiles is consistent within the error 
bars with both the self-similar scaling (S oc T) as well as 
with the observed trend, S oc T ' 65 . In any case, we would 
like to stress that we do not expect to find the self-similar 
scaling albeit we are neglecting non-gravitational processes. 




-0.2 

-0 4 




BM 




.iMr', PBM . 



-'1 11 

■ T 

i 


1 MM , , , 

1: 1 NFW 


f 


i] 


: j !ii 


k 1 MQGS 


a 

_ 


ii 






in 














hi 




ii 


PBM 






r / 9, r / R soo 

Figure 17. Same as Figure HH for the entropy profile. Entropy is 
defined as S = kTn e 2 ^ 3 and rescaled by the central values of the 
ICM density and temperature inferred from the fit to the X-ray 
surface brightness. 



The central entropy in any model scales as Topg 2 ^ 3 , where 
po oc pi and To oc p\r\. Recovering a normalisation exactly 
proportional to the central temperature would require that 
the characteristic CDM density was independent on halo 
mass (i.e. constant concentration). 

As in the baryon fraction, we argue that the character- 
istic density and radius can be more physically meaningful 
than the mass and radius at a given overdensity. Depending 
on concentration, the entropy at some fraction of -R200 can 
vary as much as a factor of 3. This can have important con- 
sequences on the S — T relation, since we have shown that 
merging systems are systematically much less concentrated 
than relaxed haloes of the same mass. A lower concentration 
implies a lower central gas density, but also a lower tempera- 
ture and a smaller value of ifeoo/n- Therefore, the net effect 
depends on the details of the mass-concentration relation. 

The S — T relation for our sample of numerical clusters 
is shown in Figure H81 The entropy, S = kT[p g /(p,m p )]~ 2/:i , 
is evaluated at 0.1 -R200 and plotted as a function of the X- 
ray temperature of the object (see Table . Systems that 
have been catalogued as mergers according to our substruc- 
ture criterium display higher entropies than relaxed clusters. 
Since mergers are more common on group scales, the S — T 
relation becomes considerably flattened when these systems 
are taken into account. 



6 CONCLUSIONS 

In this paper, we have considered four self-consistent ana- 
lytical models of galaxy clusters, based on the hypotheses 
that the hot ICM gas is in hydrostatic equilibrium with the 
dark matter halo and that it follows a polytropic equation of 
state. Two of our models assume NFW and MQGSL formu- 
lae to describe the CDM density profile, whereas the other 
two assume a /3-model for the gas distribution. One is an 



14 Y. Ascasibar et al. 



Figure 18. Gas entropy at 0.1 -R200 versus emission-weighted 
temperature. Solid squares represent relaxed clusters, empty 
squares are used for minor mergers and stars for major merg- 
ing systems. Solid line depicts the self-similar scaling, S oc T, 
while dashed line shows the observed behaviour, S oc T - 65 . 



isothermal version with /3 = 2/3 (BM) and the other is a 
polytropic model with 7 = 1.18 and j3 = 1 (PBM). 

We performed a set of high-resolution adiabatic gasdy- 
namical simulations to assess the accuracy of the basic ap- 
proximations (i.e. hydrostatic equilibrium and a polytropic 
relation). Then, we compared the radial distributions of gas 
and dark matter expected in each analytical model with the 
results for a sample of 15 numerical groups and clusters be- 
tween 1 and 3 keV. 

Our main conclusions can be summarised as follows: 



parameters. Purely adiabatic physics cannot account for a 
large isothermal core. 

(vi) /3-models fail to reproduce the gas density profile. 
While a fit to the inner part yields /3 ~ 2/3, larger values 
are required at large radii. NFW and MQGSL models are 
able to predict the average distribution, but variations in the 
polytropic index lead to a very large scatter. The asymptotic 
slope of the gas density profile is too steep compared to the 
simulation data. 

(vii) Models based on a 'universal' CDM profile are able 
to estimate the ICM properties within 30—40 per cent errors, 
when applied to the simulated X-ray surface brightness. As- 
suming a /3- model yields similar estimates for r ^ O.LR20O) 
but the shape of the inferred profiles at smaller radii can be 
severely misleading. 

(viii) Contrary to conventional SPH estimates, the en- 
tropy distribution is not a pure power law. All our analyt- 
ical cluster models predict an increasing slope with radius, 
consistent with our numerical data. The entropy profile is 
entirely determined by the characteristic density and radius 
of the CDM halo, and the normalisation depends on the 
mass-concentration relation. In particular, it does not nec- 
essarily scale linearly with gas temperature. 

In general terms, we claim that the radial structure of 
galaxy groups and clusters can be understood in terms of 
hydrostatic equilibrium and a polytropic equation of state 
for the gas, at least when only gravitational processes and 
adiabatic gasdynamics are considered. 

Although /3-models can sometimes provide a fair estimate 
of several cluster properties, they fail to provide an over- 
all description that matches all the numerical data. Models 
based on a 'universal' CDM density profile have a similar 
or better accuracy, and they are far more reliable in a qual- 
itative sense (i.e. concerning the shape of the profiles). We 
argue that these models are preferable in order to analyse 
and interpret observational data. 



(i) We find additional evidence of non-physical entropy 
losses in the standard implementation of the SPH algorithm. 
This effect can be of critical importance in the computation 
of the temperat ure and entropy profiles. The new formula- 
tion proposed bv lSpringel fc HernauislJ J2002T) appears to be 
a promising alternative to overcome this problem. 

(ii) Thermally-supported hydrostatic equilibrium is a 
valid approximation for objects classified as 'relaxed' on dy- 
namical grounds, up to ~ 0.8i?2oo- We find evidence of non- 
thermal support in merging systems. 

(iii) A polytropic equation of state with 7 ~ 1.18 pro- 
vides a fairly good description of the ICM gas. Substructure 
induces deviations from this relation. 

(iv) The density profile of our clusters can be well fitted 
by either NFW or MQGSL formulae. /3-models predict con- 
stant dark matter density near the centre, in conflict with 
the numerical results. Major mergers feature a power-law 
density profile with slope — 2 ^ a ^ — 1. 

(v) Gas temperature declines by a factor of 2 — 3 from 
the centre to the virial radius. Our projected X-ray temper- 
at ure profiles can be f itted by the 'universal' form proposed 
bv lLoken et al.l i2002f) . derived from Eulerian gasdynamical 
simulations, although we obtain different values for the free 
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APPENDIX A: CLUSTER ATLAS 

In the following image, we display our 15 simulated clusters, 
projected along the three major axes. Colours indicate the 
dark matter surface density along the line of sight. Contour 
plots represent the X-ray emission in conventional units. 

Quite interestingly, it can be appreciated that many im- 
ages of merging systems (both minor and major) appear to 
be spherically symmetric due to projection effects. Usually, a 
lot more substructure is present in the CDM density (which 
roughly corresponds to the optical emission) than in the hot 
ICM gas (X-ray). 

Indeed, we find that the gaseous component of infalling 
galaxies is very efficiently stripped, while dark matter haloes 
are able to survive during many orbits around the cluster 
centre 3 . We expect these haloes to retain some of their orig- 
inal gas when cooling is taken into account, probably under- 
going an intense star formation burst. However, all gas in a 
warm- hot reservoir is lost during the first orbit, quenching 
any subsequent star formation activity. Numerical exper- 
iments that include self-consistent cooling, star formation 
and feedback are needed in order to make a quantitative 
assessment of this scenario. 



3 A visual impression of this phenomenon can be ob- 
tained from a series of animations that can be accessed at 
http : //pollux . ft . uam . es/gustavo/VIDE0S/LCDM80/clustervideo . html 



